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Abstract 

We consider the interplay between superconducting coupling and dopant 
impurity scattering of charge carriers in planar square lattice systems and 
examine the physical conditions (doping level, temperature, local symme- 
try of coupling and scattering potentials) necessary in this model system 
to obtain a d-wave superconducting order, like that observed in real doped 
cuprate HTSC materials. Using the Lifshitz model for the disorder intro- 
duced into system by dopants, we analyze also the non-uniform structure 
of such d-wave parameter, including both its magnitude and phase vari- 
ation. The results indicate that d-wave superconductivity turns possible 
in a doped metal until it can be destroyed at too high doping levels. 

PACS: 71.55.-i, 74.20.-z, 74.20.Fg, 74.62.Dh, 74.72.-h 



1 Introduction 

The studies of the effect of impurities and defects on superconducting (SC) 
properties of metals (including SC alloys) began practically as early as the BCS 
theory had been built. Inparticular, the classical papers by Anderson |jl| and by 
Abrikosov and Gor'kov ^ indicated a substantial difference between magnetic 
and non-magnetic impurities in superconductors. If an addition of non-magnetic 
impurities has practically no effect on the value of transition critical temperature 
Tc, the presence of spin at impurity atom (leading to the Kondo effect in a normal 
metal) results in pair-breaking, that is, it transforms a singlet Cooper pair into 
an unstable triplet and rapidly suppresses Tc. In all the known cases, the s-type, 
or isotropic, SC order (apart from the heavy-fermion systems, where it pertains 

'Author to whom correspondence should be addressed. E-mail: vloktev@bitp.kiev.ua 
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to the p-type and Tc is extremely low) and, correspondingly, the isotropic gap 
near the Fermi level were considered. 

The discovery of high-Tc superconductivity (HTSC) in copper oxides faced 
the physicists to a number of problems which still remain a challenge for the 
theory. Undoubtedly, this includes the issue of HTSC mechanism, the strong 
dependencies of many (both SC and normal) properties of copper oxides on 
charge carriers concentration, specifics of weakly doped systems (first of all, 
existence of a pseudogap at temperatures above Tc), formation of stripe struc- 
tures, etc. (see, e.g., the review articles [il,|@|,||l|,|§l). Such a problem is also 
presented by the impurity effect on SC properties of HTSC systems. These differ 
from "old" or "conventional" superconductors not only in higher Tc and c?-wave 
anisotropy of the order parameter but also in the fact that here magnetic and 
non-magnetic impurities change their roles: the former are weak suppressors for 
Tc 0, Ig], while the latter (in particular, Zn substituting Cu in cuprate layers) 
lead HTSC to a fast decay 0, [|lo|, Many aspects of impurity effects in 

superconductors with anisotropic (including d-wave) pairing were already theo- 
retically studied in Refs. ||||,||l||,Q,||l||,||6l,||l3, using a range of models and 
approximations. However, these (and many other) papers did not include one 
of the most essential features of HTSC systems, the fact that they cannot be 
non-impure. 

In other words, most of theoretic approaches to HTSC are based on the con- 
cept, formulated and applied in the pioneering works one starts from an 
ideal (2D or quasi-2D) metal with given Fermi energy ep, defined by the den- 
sity of free carriers, and then consider perturbation of independently existing 
SC condensate by some extrinsic (magnetic or non-magnetic) impurities. Their 
action, local or global, affects the pre-formed and condensed singlet pairs. Of 
course, this formulation is reasonable but it lacks an essential moment for the 
conductance in copper oxides: all the HTSC's are doped metals, where (like 
the doped semiconductors) each carrier is provided by insertion of a donor or 
acceptor into the system. In its turn, this impHes that HTSC's are intrinsi- 
cally impure systems with an inherent disorder,^ and the number of impurity 
("foreign" included) atoms can not be less than at least the number of charge 
carriers. If one has in normal metals the condition k^i ^ 1 (fcp being the Fermi 
momentum, and £ the carrier mean free path between collisions with impurity 
atoms) iQ, it turns to k-pi ~ 1 in HTSC's, and they pertain to the family of 
"bad" metals with both k-p and £ defined by the doping. 

Perhaps, the first attempt to consider in a self-consistent way the charac- 
teristic tendencies for HTSC, pairing of the carriers and their localization on 
impurity atoms, was made in the authors' work It discussed the phase 

diagram of doped 2D metal in presence of s-wave pairing and showed that SC 
is not possible neither at low impurity concentration c < cq ~ e^/W (when 
all the carriers are locaHzed near impurities with locaHzation energies eo much 
less of the bandwidth W , so that cq is typically few percents) nor at too high 

^Here we don't consider possible formation of stripe structures, where ordered or disordered 
distribution of dopants can not be yet confirmed by any reliable data. 
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c (when the pair inverse Hfetime times h exceeds the SC gap). There, in gen- 
eral, the self-consistency is related either to the SC order parameter (like the 
common Bardeen-Cooper-Schrieffer or Bogolyubov-de Gennes treatments) and 
to the chemical potential. 

The present work is aimed to extend the approach by Ref. to the case 
of d-wave SC coupling and to trace the formation of the corresponding order 
parameter. It is motivated, not to the least degree, by an apparent controversy 
between the experimental evidence for d-symmetry of the order parameter in 
HTSC lH^, and theoretical claim that anisotropic pairing should not 

survive in presence of chaotically distributed isotropic scatterers [g^. For the 
sake of simplicity, we restrict the consideration to the doping range c > cq, 
where the self-consistency is only relevant for the SC order parameter, while the 
chemical potential can be putQ /i « ep ~ 3cW/A. Then, we distinguish between 
two types of impurity effects by doping. The first, the so called homogeneous 
effects, are displayed by translationally invariant single-particle Green functions 
(SPGF). They were studied earlier by various means |Q but, as a rule, intro- 
ducing the disorder through a single parameter Va of Anderson's model 
In contrary, we employ the Lifshitz' model of disorder |p5|| , characterized by 
two independent parameters: c and the impurity potential Vl- They produce 
an equivalent Va ^ \/c{l — c)Vl, but not vice versa. Within this model, more 
adequate for doped HTSC systems, we conclude about persistence of d-wave or- 
der parameter under homogeneous impurity effects. Also we explicitly consider 
other type of effects, inhomogeneous, due to local variations of order parame- 
ter near impurity centers. This involves two-particle Green functions (TPGF), 
besides usual SPGF, and yields in possible limitation for SC at high enough 
dopant concentrations. 

At least, we would like to appreciate this great honor and also a pleasant 
opportunity for us to publish this work in the Low Temperature Physics issue 
dedicated to the memory of outstanding physicist L.V. Shubnikov, whose con- 
tribution to the low temperature physics in general and to superconductivity in 
particular can not be overestimated. 

2 Hamiltonian and Green functions 

We start from the model electronic Hamiltonian in band representation 

^^Y.{Y. ^k4,^ck.<x -jjYl [^7k7k'4,|clk,|C_k',iCk',T- (1) 

k I o- k' 



^However, it is known tha 
of very low doping (see, e.g. 



,n differ essentially from the Fermi energy ep in the limit 
)• 
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where Ck,cr and c]^ ^are the Fermi operators for a charge carrier with wavevector 
k and spin a. The simplest band energy = At — 2t (cos ak^ + cos aky), with 
full bandwidth W = 8t,is expressed through the amplitude t of carrier hopping 
between nearest neighbor sites|^ (of the total N in the lattice with a constant 
a). The parameter V models the attraction between two carriers with opposite 
spins on such sites, the factor 7k = (cosafcx — cos aky) 9 (e^ — has the d- 
wave symmetry and is effective only for quasiparticle energies = — £f 
smaller than the "Debye energy" Ed- The latter is understood as a characteristic 
energy of intermediate (FroeHch) boson, and in what follows we suppose the 
condition £d < to hold and a BCS-shell to exist (the alternative ^ < Ed, 
possible for underdoped HTSC systems, will be considered elsewhere). The 
impurity perturbation Vl expresses the shift of on-site electronic energy on a 
random dopant site p, where the negative sign takes an explicit account of 
carrier attraction to the ionized dopant and, for simpHcity, we consider this 
perturbation localized on a single site. With usual BCS ansatz: c_kaCk,f = 
(c_k4Ck,|) + <y5k, and in neglect of quadratic terms in pair fluctuations (pit:, Eq. 
|l| leads to a bilinear form for H' = H ~ /iiV: 

CkCk ,^Ck,o- - (AkC-kaCk.T + /^•c•) - (2) 
E e-'('''-''>^et,,Ck,. . 

p,k',cr 

Here the gap function is deflned by the self-consistency relation 

Ak = (ck',TC-k',i) , (3) 

k' 

extending the common BCS gap equation to the d-wave case. A non-uniform 
system can be treated within the formalism used formerly for impurity problems 

in SC passing to the Nambu spinors i/'k = (cj^ pC^k.ij and i/'k, and 

deflning the Fourier transformed matrix Green function (GF) 

Gk,k' (e) ^« V'klV'l. »e= r e^(^-'0)*({^k(0,4'})'^^- (4) 

Here A denotes a 2 x 2 matrix in Nambu indices, (...) is the quantum sta- 
tistical average, and {a (t) , b (0)} the anticommutator of Heisenberg operators. 
In the GF's below we omit their explicit dependence on energy e but distin- 
guish between their diagonal and non-diagonal forms in Nambu (N) and mo- 
mentum (M) indices. Then, applying the Heisenberg equation of motion (EM) 
indi}k/dt = [H^tpk] in Eq. |, we arrive at the Dyson's type EM for SPGF's: 

Gk,k' = 4,k' - y "'^""""^'^Gk^k' (5) 

P,k" 

■'So, we do not take into account next neighbor hoppings. 



k 
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where the non-perturbed SPGF g'^'' = (e — ^k^^s — + iO) ^ , and the scat- 
tering matrix V = Vj^t^ include the PauU matrices Ti. 

For a disordered system, the relevant (observable) characteristics are de- 
scribed by the so-called self-averaging GF's, whose values for all particular real- 
izations of disorder turn practically non-random, equal to those averaged over 
disorder ||2^. The most important example of such a function is the M-diagonal 
SPGF, Gk = Gk,k-The general solution for Eq. ^ in this case can be written 
(see Appendix E) as 

^ 1 



G 



(0) 



-Si. 



(6) 



where the self-energy matrix Sk is given by the so-called fully renormalized 
group expansion (GE) li|,|l4l,|§ 



Ek = -cV 



1 + GV 



n=iO 



-ik-n 



+ ^On^nO 



1 — ^On^nO 



(7) 



Here the integrated SPGF matrix G — N ^ Sk^k, and the matrices Aq-^ of 
indirect interaction between scatterers at sites and n are 



k'/k 



^Gk' 



l + GV 



(8) 



The rectriction to k' k at single summation in Apn should be complemented 
by k" ^ k, k' at double summation in the product ^on^no, but such restrictions 
can be already ignored in ^on^on^no and higher degree terms ||2^, resulting 
from expansion of the r.h.s. of Eq. 0. 

Many observable characteristics of SC state follow from GF's, using the 
spectral theorem representation 



{ab) = 



de 



-Im < b\a >e 



(9) 



3 e^'^^-^) + 1 

where the chemical potential ^ is defined by the overall electronic concentration 



-E 

k 



de 



e/3(e-M) + 1 



ImTrTsGk- 



(10) 



On the other hand, c is just the concentration of dopant centers: c ~ 1, 
which give rise to carrier scattering, and the carrier concentration only gets close 
to (but never exceeds) c in the regime of doped metal, for c above a certain met- 
allization threshold cq (for quasi-2D dispersed £k, it is cq ~ exp {—-kWIAVl) ^ 1 
[E9)). At this condition, the self-consistency implied by Eq. 10 is not necessary. 
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and a good approximation^ for chemical potential is fj, ~ 3cW/ 4 (see Appendix 
A). Then the gap equation, Eq. ^, presents as 

Ak = -Sr y r -^r^^ 7lm TrfiGk' , (11) 



2N ^ 7-00 c/'(^-A') + 1 

and its solution discussed in Appendix C for the uniform case (Vl = 0) is simply 
Ak = A7k, with the ratio r = 2A/A;b7c being e^/^ times the s-wave BCS value 
TBCS ~ 3.52. 

Another important self-averaging quantity is the integrated SPGF matrix G 
itself, since the density of states p (e) is just 

p(e) ilmTrG. (12) 

TT 

For a non-perturbed system, Vl ^ 0, G ^ G*^"^ = N^^ g'"^\ and calculation 
of the imaginary part of G*^°-' within the nodal point approximation (Appendix 
B) leads to the standard d-wave density of states: 



p(e) p(o) (e) = ilniTrG(°) = ^arcsin min ( 1, 



(13) 



where po ~ 4:/{nW) is the normal Fermi density of states of doped (quasi-2D) 
metal. Respectively, the real part of G'^Hs 



ReG(°) = spo —r^ — r - ^61 (A - e) arccosh— 
Then Eqs. |l^ and ^ can be unified into a single analytic form: 



(14) 



W TT / A TT 

arccosh i 



p{W-p) A\ e ^2)y ^^^^ 

since at £ > A one has arccosh(A/e) = z[7r/2 — arcsin(A/e)], thus restoring 
Eq. |l^. But it is just the growth of (real) arccosh term at e < A that permits 
existence of a low-energy (Sres ^ A) resonance feature in Re(l + G^)^^, and 
hence in p(e). Such a resonance was formerly discussed for a d-wave SC with 
low enough concentration c (so that G w G^"^) of "foreign" impurities producing 
strong enough perturbation Vl |[14| , and it is similar to the known low-frequency 
resonance by heavy impurities in acoustic phonon spectra However, in the 
situation of our interest here, when both Vl and c are not small, G can be 
essentially modified compared to G*^°' and this is expressed in a very complicate 
way by Eq. 0. To simplify the task, certain self-consistent procedures, like the 
CPA method, quite useful in the theory of normal metals can be employed. 
Similar approach was formerly proposed for an s-wave SC doped system ||l9|, 
and here we begin with the analysis of a self-consistent solution for Eq. |H| in the 
d-wave case. 

*This approximation is actually justified by the fact that for c > cq the Fermi level ep 
of metallic phase is well higher than the conduction band edge, and one can hardly suppose 
existence of local pairs and related inequality < ep at these concentrations. Therefore, in 
what follows we do not distinguish between and £p. 
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3 Uniform doping effects in self-consistent ap- 
proach 

If the GE series, Eq. ^ is restricted to its first term, the self-energy matrix Sk 
turns in fact independent of k: 



S ^ -cV 



l + GV 



(16) 



and substitution of Eq. into Eq. ^ defines the self-consistent approximation 
G^'"^ for M-diagonal SPGF: 



Gi 



G 



(0) 



g(sc) 



S^^'^) = -cV 



[l + g(" 



(sc) 



(17) 

(18) 
(19) 



To solve this system, we first parametrize the self-energy matrix, Eq. |l^: 

= Eo + + S3T3, (20) 

Ei being some complex-valued functions of energy. Then integration in Eq.p^ 
within the nodal point approximation (Appendix D) results in 



G(^=) = Go + Gin + Gafa 



with the coefficients 



(21) 



Go = (e - Eo) po 



— - arccosh — h arccosn — nr 

2A V £ - So £ - So 



^i{W-^l) 



(22) 



Gi = Eipo 



(e - I]o)2 - (A + Ei)^ + ./(£ - Eo)2 - (A - 



W 



G3 = Po 



^(H^-/i) 



(23) 



(24) 
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Substituting Eq. ^ into Eq. |l^, we arrive at 

cli [Fl (Go + Gin) - (1 + VlGs) fg] 



(1 + VlG3 



(25) 



25 



and g2|, 
or that 



!J with Eq. |2^, we immediately conclude 
is in fact N- diagonal, which is extremely 



Comparing Eqs. 
that El = Gi = 0, 
important. Physically, this means that (within the self-consistent, Hnear in c 
approximation) the scattering by dopants does not influence the d-wave order 
parameter, and this can be directly related to the fact that the s-symmetry of 
impurity perturbation Vl is orthogonal to the rf-symmetry of SC pairing V. It 
also applies to more realistic models of dopant perturbation in HTSC (e.g. with 
plaquette- or dumbbell-like anisotropy fsj), as far as their symmetries do not 
coincide with that of the order parameter. CompHcations arise when they do 
coincide, as was found for isotropic perturbation on s-wave order with all three 
Si being non-zero jl9|, hence the apparently "harder" d-wave system turns in 
fact "easier"! 

Using the fact that Si 



and the relation cosh(a; 
is brought to a very simple form: 



■i7r/2) 



Eq. 



A . f Go 

= sm a — 



e-Sn 



A 



npo e - 



(26) 



with a = VFA/[7r/i (W — ^)] <C 1, while the same comparison for the two non- 
zero components of S^'^'^^: Sq and S3, gives: 



So 



(1 + VlGs)^ - Vj^G'o 



cVlGo, 



(1 + Vi^GsY - VlGl = -cVi^ (1 + FlG 



(27) 



(28) 



From Eq. ^ we estimate IS3I ^ cVl, hence within the relevant energy region 
|e| <C Ed, the function G3 from Eq. U is reasonably well approximated by a 
(negative) constant = pa \ii[p/{W — p)]. Then Eq. |2^ turns quadratic for Go- 



Gq -I- -^Go 
^0 



= 0, 



(29) 



with V = Vl/(1 + Vi^gs) w 14 ln(l/co)/ ln(3c/4co). The system, Eqs. ||, ||, 
fully defines the self-energy Sq and other uniform physical properties of a dis- 
ordered d-wave system, and its solution can be found (in principle numerically) 
within the whole relevant energy range. It turns especially simple if |So| ^ cV 
(this proves to hold at least for e ^ Sres), then the proper solution to Eq. ^ 
is: Go ~ So/cT^^, and from Eq. ^ we obtain the following equation for a single 
important function So: 



So 



sm a 



/3So 
e-So 



(30) 
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with P — A/[ttcV'^pq] ^ 1. It defines the self-consistent density of states 

pise) ^ iiniTrG(-) (s) = (31) 

at lowest energies. This approach results free of infrared logarithmic divergen- 
cies, appearing in the integrals of perturbation theory and thus permits to 
avoid applying heavy field theory methods for white-noise scattering potential 
[ p^ , whose adequacy to the case of discrete random dopants is not clear. 

The exact value of the density of states at the very center of the gap, p (0) = 
p{e ^ 0), is also of a particular interest in view of the known claim about 
existence of a non-zero "universal" value p (0) ~ c/ln(l/c) if Vl is sufficiently 
strong iQ, However, we conclude from Eq. ^ that in the limit e ^ 0: 

^ ln(2A//3|e|) 
and hence the self-consistent density of states at e ^ 0: 

(s) ^ (321 

vanishes even faster then the non-perturbed function p'"-* (e), Eq. This 
produces a certain narrow "quasi-gap" (not to be confused with the pseudo- 
gap observed at T > Tc in the underdoped regime) around the center. For 
comparison, the estimated p (e) from the two first terms of Eq. tends to zero 
Hnearly in e with corrections ~ [Q, while the field-theoretical analysis |Q 
predicts p (e) e'^ , with the non-universal exponent being (in our notations) 
CO = tanhln y/ tt^AW/2cV^, that is always < 1 and even not excluding < 0. 

The discrepancy between our results and the before mentioned "universal" 
behavior originates in the improper use of unitary limit Vl ^ 00 at neglect 
of the 1 -I- VlGs term in Eq. leading to the relation So — —c/Gq. But 
the true limiting relation is inverse: Eq = cV'^Gq, with V finite for Vl 00, 
and also the unitary limit fails at any^ finite Vl when e ^ 0. Finally, the 
existing experimental data do not confirm the "universal" p (0) value, but seem 
to favor the conclusion about existence of a strong low-energy resonance in 
p{e) with possible quasi-gap at the very center ||3^, though experimental 
observations at such low energies of course need extremely low temperatures. 

Note however that the self-consistent treatment of the low-energy resonance, 
at Eres Aln(3c/4co)/(7rln27r) for the case of "own" impurities (Fig. ??), 
requires already solution of the full system, Eqs. and, in view of a 

probable underestimate of this hump (like that in normal systems [|8|,[Q), it 
should be better described by the exact GE, Eq. ^. 

The obtained E^'"^) can be in principle directly inserted in Eq . p7| , in order 
to use the resulting G^**^' for correction of the gap equation, Eq. However, 
at important there quasiparticle energies ~ ed, renormalization effects are 
negligible, and thus A remains well approximated by the result of Appendix C. 

^Unless the special case l/V = corresponding to c ^ 4co/3, while the actual consideration 
is for c ^ CO . 
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4 Non-uniform effects 



The SPGF's, considered in the previous section, describe the uniform self- 
averaging characteristics of SC state. The next important question is the be- 
havior of fluctuations of the order parameter (both its amplitude and phase) 
in an inhomogeneous system, which should be closely related to the breakdown 
of superconductivity in the overdoped regime. A strong local suppression of 
d-wave order close to a single "foreign" impurity was predicted theoretically 
and observed experimentally ||ll| . In the general case of finite concentration of 
scatterers, the local d-wave order can be characterized by the operator 

= ^ E 7ke<'^-''') "c_k,iCk',T, (33) 

k,k' 

such that its mean value (generally coniplex) defines the uniform gap parameter: 
N-^ En (^n) = N-^ Ek7kAk = A Q. On the other hand, it is natural to 
characterize local fiuctuations of order parameter by the mean squared disper- 
sion of fin (identified with the dispersion of gap parameter): 



6' 



N 

Using Eq. we write 



^E((^n)-(^^n)^). (34) 



= ^E E e^^''^+''="''^"''^^-"7k,7k. [(c_kuCk3,TC_k,aCk„T)- 

n ki,k2,k3,k4 

(35) 

- (c_ki,iCk3,T> (c-k2,iCk4,T)] > 

or, summing up in n to close the sum of 4 momenta, and using the spectral 
representation, Eq. 0, for T = 0: 



ki,k2,q 



r poo 

E Tki7k2 / delm < c_ki,iCk2,T|c-k2+q,iCki-q,T > - (36) 

1 ,k2,q 

/>00 />CXD 

- / delm < c_kia|ck2,T > / c?elm < c_k2+q,i |cki-q,T > 
Jo Jo 

Here, besides the before used SPGF's, the TPGF ^ c_ki,j.Ck2,T|c-k2+q,iCki-q,T ^ 
appears. An explicit calculation of this function for non-perturbed (Vl = 0) SC 
system (see Appendix E) gives the following result 

« c_k„iCk2,T|c-k2+q,iCk,-q.T »= So,^ El^^^^Elf^- 'ill El • ^^^^ 

Here E^ = Ck + ^kj ^ind three other terms in numerator are not mutually odd in 
7ki and 7k2 and thus do not contribute into S^. Then, it is easy to see that the 



10 



resulting contribution into Eq. |3^ from M-diagonal (q = 0) and N-non-diagonal 
TPGF's by Eq. |3^ exactly cancels with that from SPGF's, which confirms 
uniformity of d-wave order in this case. The most important contribution to 5"^ 
at 14 7^ comes from two consecutive scattering processes in the l.h.s. of an 
N- and M-non-diagonal (q 7^ 0) TPGF, first c_ki,j,Ck2,T c_ki+q,j,Ck2,T and 
then c_ki+q.iCk2,t c_ki+q,xCk2-q.T (or vice versa, Fig. ??), on the same 
scattering center p. It is linear in c, while the contribution from SPGF's in 
this order, —2/\V Jq°° delm TrriS, is zero, accordingly to N-diagonal form of E 
estabHshed in previous Section. Hence we generally estimate the gap dispersion 
to grow with c as |(5| = Ay^c/ci, where ci ^ (V^/Vl)^ defines the upper critical 
concentration for SC at T = 0. 

A more detailed analysis, resolving amplitude and phase fiuctuations, can 
be done in a similar way, but considering separately two operators: 

f^n,+ = — and VL^ - = , 

such that their mean values (f^n,±) lead to real and imaginary parts of the order 
parameter, and constructing the corresponding dispersions: 

.2 1 



This approach should be particularly important at extension of the theory to 
finite temperatures, in order to establish the dominating type of fiuctuations 
due to the static disorder, responcible for breakdown of SC order at T = Tc, 
and its possible role in the persistence of a pseudogap in the density of states p 
at T > Tc. 



5 Conclusions 

The above presented analysis shows that the disordered structure of doped 
HTSC systems is crucial for many of their characteristic properties and for 
existence of SC order itself. The interplay between doping and disorder effects 
can be briefiy resumed as follows. Superconductivity onsets with the metalliza- 
tion of the system, at critical concentration c ~ cq resulted from the competition 
between kinetic energy of charge carriers in regular lattice and their attraction 
to random dopant centers. The uniform d-wave order parameter A develops 
with growing number of charge carriers as A '--^ -Jc |^ and saturates at certain 
optimum doping Copt ~ ed/W^, when the relation /j, > Ed becomes to hold. With 
further growing c, the increasing local fiuctuations of A bring it to collapse at 
some upper critical concentration ci ~ (V^/Vl)^, resulted from the competi- 
tion between pairing and scattering potentials. This picture is quantitatively 
satisfied with a very natural choice of parameters W ^ 2 eV, Vl ^ 0.5 eV, 
V ~ 0.22 eV, Ed ~ 0.2 eV, giving plausible estimates: cq ~ 5%, Copt ~ 15%, 
Tc,max ~ lOOK, ci ^ 20%. The forthcoming work should also specify such im- 
portant moments, left beyond the scope of this consideration, as the disorder 
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effects on the cusp of density of states p{e) at e = A, the matching conditions 
between the self-consistent and GE descriptions of SPGF, the exact numerical 
coefRcient for the critical value c\, etc. And, of course, it is of principal interest 
to extend the present self-consistent treatment to the case of finite temperatures 
up to Tc, in order to obtain a quantitative estimate for the bell-hke Tc(c) shape, 
and further to T > Tc, to study the role of doping disorder vs rf-wave SC cou- 
pling in the formation and subsequent merging (at c ^ Copt) of the pseudogap 
in normal density of states. 

This work was supported in parts by Portuguese Program PRAXIS XXI 
under the Grant 2/2.1/FIS/302/94, by NATO Program OUTREACH under 
the Grant CP/UN/19/C/2000/PO, and by Swiss National Science Foundation 
(SCOPES Project 7UKPJ062150.00/1). One of us (V.M.L.) expresses his grat- 
itude for a warm hospitality he felt at visiting Centro de Fisica do Porto. 



6 Appendix A 

For the normal phase, we consider the dispersion law £k = 4:t—2t (cos ak^ + cos aky) 
and relate the dopant concentration c to the number (per unit cell) of occupied 
states below /u: 



2^2 r 2ct^ /'arccos(l— /■arccos(2— 4/^/1^— cos a/cg;) 

^ / dk = — ^ / dkx / dky 

{2nf Je.<^ Jo Jo ^ 



(38) 



du f dv p ( 



Here the dependence of the integral 

'■^ arccos (2 — x — u)du 



F{x) 



l-x 



is very well approximated by a simple linear function F (x) ~ 5x/'S in the 
whole physically important range < x = 4fj,/W < Xmax (where Xmax ~ 0.6 
corresponds to fimax ~ 0.15VK at maximum physical doping Cmax ~ 0.2). Then 
we readily arrive at the estimate 

37r2 3cW 
'^ 40 4 ' 

cited above. 



7 Appendix B 

For the uniform SC system with d-wave gap, we perform integration over the 
Brillouin zone with the parametrization: k— k, = a~^{po^ei+A~^rieiXez), e, = 
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(±l/\/2, ±l/\/2, O), = (0, 0, 1), near 4 nodal points — arccos (1 — iijW) 
of the gap function Ak — ^1^^ ■ This integration for G^'^Hurns into 



N ^ 



k 

(0) _ 



k 

— ED 



W-fj. 



(39) 



where we have defined the matrix function 



■qTi 



^2 _ ^2 



The integration in ^ (normal to the Fermi surface) in Eq. ^treats the BCS shell, 
[— eDj^d], separately from the out-of-shell segments, [—^,—£0] and [eD,W^ — m], 
where the gap parameter A turns zero (together with 7k) and no integration in 
rj — A7k is needed. Eq^||9| permits to define explicitly the coefficient functions 
gi in the general form G'-">= go + giTi + (73X3. Let us denote = — 77^, then 
the shell contribution to go result from the integral: 

r'^ 2 f z \ 2 
— = — arctann V vn \ k, — 



2iiT 

z 



which is followed for its last term by: 

dr\ 



A 



A 

= — 2arccos — . 
2 e 



The out-of-shell contributions are: 

If , e , £ 

— = — arctanh arctanh — 

^ £^ - £ V M 

and 

d^ 1 



1 
At 



1 

Ed 



£2-^2 W-fi 



1 



To find gi and gs, we use the evident equalities 

rjdri = 0, 



and: 



(40) 



(41) 

(42) 
(43) 



2 _ ^2 



Summing up Eqs. ^0[Q, we obtain 

W 



go = £po 



1 ^2_e2 

7: in T 

2 (W- - £2 



In- 



W - n 



/XX. ^ — -r i arccosh — 
^i{W - fi) A V e 



ITT 

T 



in accordance with Eq. IT^, gi = 0, and 513 = po ln[^/ {W — fi)]. 



(44) 
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8 Appendix C 



The gap equation, Eq. ^ for uniform (Vl = 0) rf-wave system at T = 
transforms into: 

1 16 f°° de 

-Im / 7fdfj / / 5 5 (45) 



A 7rA3 Jo Jo Jo e^-e-V^~^0 

with A = Vpo. The sought quantity is the gap ampHtude A. On the r.h.s. of Eq. 
^we perform the elementary integration in e, using the relation Im {x — iO)~^ = 
ttS (x): 

16, /"^ o , r° r de 



A 



8 / ffdri 







/o VFT^ 

Then, integrating out in ^ and passing from r] to y — huju/rj, we present Eq. Eel 



as 



arcsinhy dy 8 Ed Sd vEt 
= -arcsmh— + 4 — - 



£D y4 3 A 3A2 



4 /eD\3 . , A 8 / 2£D 1 

— — arcsmn — sa — In — : — 

3VA/ £d3VA 3 

The equation for critical temperature Tc, corresponding to A = 0, reads in this 
case 



1 16, 
— = — Im 
A 3tt 



Jo \2kBTj Jo e^-e-^^ 3 \TrkBTcJ 



with the Euler constant 7e ~ 1.781. Hence for the d-wave case the effective 
coupling constant is A — 8 A/3, that is 8/3 times the "Hamiltonian" value A, 
which can serve as one more explanation for high Tc itself. Also, the ratio 
r = 2A/kBTc results here exp (1/3) times the common s-wave BCS value = 
27r/7E « 3.52, reaching as high value as Vd « 4.92. In its turn, this means 
that, for equal other conditions (say, po and V), s-condensate turns out more 
stable to thermal fluctuations and requires a higher Tc to destroy it than d- 
condensate. Of course, this is directly related to the absence of gap in the latter 
case, permitting quasi-particles to exist at any T < Tc. 

9 Appendix D 

Calculation of the self-consistent SPGF G*^*^-' — Go + GiTi + G3T3 generalizes 
the scheme of Appendix B: 



14 



N ^ 



k - Po 



(47) 



de?(e-Eo,e + S3,Si) 



de?(e-So,e + S3,Si) 



Next we set — [e — Sq)^ ^ + ^i)^ and z = r + ^7, and pass from ^ to 
a; = ^ + S3, where S3 can be taken real (as seen e.g. from the final result Eq. 
H). Then the analogue to Eq. |^ is: 



dx 



£1 _ j,2 2z 




(£d + S3 + rf + 72 



(£D - S3 - rf + 72 



(eD-S3 + r)V7' 



£D + S3 + r 
arctan h arctan 



(£D + S3-r)'+72 
Ed + S3 - r 



7 



Ed - S3 + r eo - S3 - r 
+arctan h arctan 



^ _ . 7r + 27/£D 
Ed z 



(48) 



where the small term 27/eD can be safely dropped. The next integration, in 
2/ = ?7 + Si, is done only on in/z term accordingly to: 



A+Si 



dy 



A + Si A -Si 

arccos h arccos- 



(49) 



/-A + Si V(e-So)2-2/2 £ - So £ - So 

which is relevant for Go, supplemented with 

r^+^i ydy _ 

4ASi , , 

; ^= (50) 

'(£ - So)2 - (A + Si)' + ^(£ - So)2 - (A - Si)' 

for Gi. The out-of-shell integration of the mentioned components is much easier, 
giving: 



dx 



1 



£-Sn 



arctanhi^-arctanh-^-^" 



M - S3 



£D — S3 
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and 



1 

dx 



1 

1 



1 



(51) 



(52) 



/,^_S3 (£-So)2-a;2 M/-Ai-S3 £0-^3 

Here also E3 can be disregarded beside VF,/i,eD, then the two — I/ed terms 
cancel out with that from Eq. 48, Now, combining Eqs. Eq-p2l we obtain 



Go = (e - So)po 



W 



— arccosn — h arccosn — 



and 



2A 



W 



At least, G3 is obtained after 



(A + Ei)^ + V(e 



1 (£D - S3 



= - In 



D+S3 



1 

2A 



xdx 
— x"^ 

xdx 

A+S 



2 (eD + S3)'-z2' 

= 1 (ep + S3)" - 
2 (VK-/i + E3)"-2 

= iln (a^-Ss)'-^' 



(£D-S3)'-z2' 



-A + Si 



A2 
3 



y2 



in the form 



G3 = po S In 



W -11 



+ 2E; 



(53) 



(A-Si)^ 



(54) 



10 Appendix E 

Search for solutions of Eq. || in the form of GE consists in consecutive itera- 
tions of its r.h.s., separating systematically the GF's already present in previous 
iterations ||2^. Let us start from M-diagonal SPGF Gk, than the iteration se- 
quence begins with singling out the scattering term with Gk itself from those 
with Gk',k, k' 7^ k: 

Gk = Gf + Gf 1 e^('^-''')■''^Gk^k = 
k',p 
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= GL"'+cG("'yGk + G("'^^ E e<''-''')■PGk^k. (55) 

k'#k,p 

Then for each Gk',k we write down again Eq. || and separate the scattering 
terms with Gk and Gk',k in its r.h.s: 

Gk^k = G('!)yl5:e<•''-•'")■PGk",k = 
k",p' 

= cG(?i>Gk',k + GL?y^e<'''-'') PGk + g'^,^V^ ^ e^'^'-'') ^ Gk+ 

p'#p 

+ ^k'^^^ E c<'''-''")p'Gk",k. (56) 

k"7tk,k':p' 

Note that the p' = p term which gives the phase factor e*('' ~'') p in the r.h.s. of 
Eq. coherent to that already figured in the last sum in Eq. |55|, is explicitly 
separated from incoherent ones, 

e<k'-k).p' 

p' ^ p (but there will be no such 
separation when doing 1st iteration of Eq. || for M-non-diagonal SPGF Gk',k 
itself). Continuing the sequence, we shall explicitly collect the terms with the 
initial function Gk, resulted from: i) all multiple scatterings on the same site p 
and ii) on the same pair of sites p and p' ^ p. Then summing of i) in p produces 
the first term of GE, and, if the pair processes are neglected, it will coincide with 
the well known CPA result |^^. The second term of GE, obtained by summing 
of ii) in p, p' ^ p, contains interaction matrices ^p',p generated by multiply 
scattered functions Gk',k, k' ^ k, etc., (including their own renormalization) . 
For instance, the iteration for a function Gk",k with k" ^ k, k' in the last term 
in Eq. ^ will give: 

Gk",k = GL°)t>l e<''"-''"')■''"Gk"^k = 
k"',p" 

= G|°^f>le<''"-'')PGk + G(°]yle<''"-'')p'Gk+ 
k TV ^ N 

+ terms with Gk'.k and Gk".k + terms with Gk"',k (k'" 7^ k, k', k"). (57) 

Consequently, the GE for Gk obtains the form given by Eq. ^. 

Now turn to TPGF <C c_ki,iCk2.T k-k2+q,iCki-q,T ^) beginning from EM's 
in absence of scattering which develop into a 4 x 4 matrix form: 

B(a,C2,Ai,A2)/- Vqrf, (58) 



B(6,6,Ai,A2) = 



/ e 


-6-6 


Ai 


A2 









Ai e- 


Ki - 6 





A2 






A2 





£ - 6 + 6 


Ai 




[ 





A2 


Ai e - 


Ki4 


-6 ) 
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with 4- vectors: 



/ = 



/ < C-ki,iCk2,T|C-k2+qaCki-q,T > \ 
< C+ .pCk2,T|C-k2+q,iCki-q,T > 
< C_ki,XC^k2,il'^-k2+q,iCki-q,T > 
V < C+ ^C+k^ Jc_k2+q,iCki-qJ > / 



/ 2A1A2 \ 

Ai 
A2 

V 2 ) 



and ^1 = Cki , ii = Ck2 , ^1 = ^ki , A2 = Ak2 - The solution to Eq. ||: / = 
^0 cfi^^ (^1,^2,^1,^2) A, gives just the result of Eq. ^ for the 1st component 
of/. 

In presence of scattering, we consider only the M-non-diagonal (q 7^ 0) case, 
then Eq. turns into: 



(59) 



where the vectors of "single scattered" TPGF's are 



I < C_ki+qaCk2,T|C-k2+q,iCki-q,T > \ 
C+ _q^Ck2,T|C-k2+q.iCki-q,T > 
< C_ki+qaC^k2,il"^-l'2+q,iCki-q,T > 



V < C+ -q,Tcik2,il^-k2+qaCki-q,T > / 

/ < C_kiaCk2-q,T|C-k2+qaCki-q,T > \ 
< C^l^|Ck2-q,T|C-k2+q.iCki-q,T > 
< C_ki,iC^k2+qal'^-k2+q,iCki-q,T > 

V < C+ ,TC^k2+q,il"=-k2+q,iCki-q,T > / 

and the 4x4 matrices: 



f'2 



Ax 



Next, EM'S for /{ 



V 
V 



A. 



1^?3 

-v% 



i?(C;,6,A;,A2)/( ^--e-'^iPAi/", 



(60) 



i?(Ci,e2,Ai,A^)/^ = --e-"iPA2/", 



N 



with £,[ = Cki-q, C2 = Ck2-q, A'l EE Akj-q, A^j = Ak2-q, coutaiu the "double 
scattered" TPGF /" =<C c_ki+q,j.Ck2-q,t|c-k2+q,i'^ki-q,T which is already 
M-diagonal and hence can be taken just in the form of Eq. Finally, the 
solution 



/ = (6,6, Ai, A2) {(5o,q + ^ [liB-1 (6, 6, A'l, A2) 
+I2B-1 (a, ^2, Ai, A'2) I2] (6, ^2, A'l, A'2)} d, 



(61) 



18 



defines the contribution ~ c{Vi^A/V)'^ into (5^, Eq. the factor being due 
to Ai,A2-odd terms from (^i, ^2, Ai, A2) and due to dominating, 

zeroth order in Ai, A2, terms from B'^ (^i, ^2, A'^, A2) B"^ (^;,C2, A'^, A^ and 
5-i(a,e^,Ai,AyB-i(el,e^,A;,A^). 
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2,0 




Figure 1: Density of states in a d-wave SC system (heavy solid line). Thin 
solid line shows the non-perturbed p^°^ (e), arrow indicates the solution for 
Re[(l + FlGs)^ — V^Gl] = corresponding to the low-energy resonance (other 
solution at e « A is not resonant) at the choice of parameters W = 2 eV, 
Vl = 0.5 eV, c = 0.1. 
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Figure 2: Schematic of consecutive scattering processes (see Appendix E), bring- 
ing an M-non-diagonal TPGF /, through intermediate, partly M-diagonal f{ 
and /j, to fully M-diagonal /": sohd vectors standing for l.h.s. operators get 
fully nested with dashed vectors for r.h.s. 
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